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ABSTRACT 

XMMUJ004243. 6+412519 is a transient X-ray source in M31, first discovered 2012 January 15. 
Different approaches to fitting the brightest follow-up observation gave luminosities 1.3-2.5 xlO 39 erg 
s _1 , making it the second ultraluminous X-ray source (ULX) in M31, with a probable black hole 
accrctor. These different models represent different scenarios for the corona: optically thick and 
compact, or optically thin and extended. We obtained Chandra ACIS and HST ACS observations 
of this object as part of our transient monitoring program, and also observed it serendipitously in a 
120 ks XMM-Newton observation. We identify an optical counterpart at J2000 position 00:42:43.70 
+41:25:18.54; its F435W (~ B band) magnitude was 25.97+0.03 in the 2012 March 7 observation, 
and >28.4 at the 4cr level during the 2012 September 7 observation, indicating a low mass donor. We 
created two alternative lightcurves, using the different corona scenarios, finding linear decay for the 
compact corona and exponential decay for the extended corona; linear decay implies a disk that is >5 
magnitudes brighter than we observed. We therefore favor the extended corona scenario, but caution 
that there is no statistical preference for this model in the X-ray spectra alone. Using two empirical 
relations between the X-ray to optical ratio and the orbital period, we estimate a period of ^9-30 hr; 
this period is consistent with that of the first ULX in M31 (18lg hr). 
Subject headings: x-rays: general — x-rays: binaries — black hole physics 
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1. INTRODUCTION 

Ultra-luminous X-ray sources (ULXs) are point sources 
with X-ray luminosities exceeding the Eddington limit 
for a 10 Mq black hole (~1.3xl0 39 erg s" 1 ) that 
are unrelated to g alaxy nuclei (see e.g. iRobertsI 120071 : 
IFeng fc Sori^lWlt for reviews). ULXs arc often asso- 
ciated with regions of high star formation such as spiral 
arms of galaxies, and therefore many are l ikely to be re- 
lated to high mass X-ray binaries (see e.g. iMapelli et al.l 
2010; Min eo et al.l 120121 and references within). Some 
ULXs may contain intermediate mass blac k holes (ESQ 
243-4 9 HLX-1 being a strong candidate, iFarrell et al.1 
2009). Other ULXs may harbor stellar mass black 
holes, exhibiting true super-Eddingto n accretion (see e.g. 
iGladstone et al . 2009; Barnard 2010), or app arent super- 
Eddi ngton luminosities due to beaming (e.g. King et all 
12001 . Since ULXs are simply defined by their observed 
luminosities, their membership is likely to include a wide 
range of diverse systems. 

High quality ULX spectra are often described by a two- 
component emission model, consisting of a disk black- 
body component, and a component such as COMPTT 
that represents inverse-Comptonization of cool photons 
on hot electrons in an accretion disk corona. Some peo- 
ple expect the temperature of the seed photons (To) to 
equal the inner disk temperature (kT; n ), assuming that 
the coro na is compact and so only sees the inner dis k 
(see e.g. IRoberts et all 120051: iGoncalves fc Sonal 120061) . 
However, many hi gh quality ULX spectr a reject models 
where T = kT; n (Gladstone et al.l [2009^1: instead, suc- 
cessful models fall into one of two types. 



The first type resu lts in cool, opti c ally t hick coro- 
nae where T > kT in ; Gla dstone et al.l (|2009D proposed 
that these spectra represent a separate "ultra-luminous 
state" where an opaque corona hides the inner disk 
from view, meaning that the measured value for kT; n 
comes from further out in the disk, outside the corona. 
ULXs fitted this way tend to have a s oft excess, and a 
roll-o ver in the spectrum above 3 keV (jGladstone et al.l 
2009). This state is distinct from the spectral states 
observed in Galactic X-ray binaries, where the optical 
depth is <1. The corona must be compact, so that we 
still see the inner disk in X-rays; this is also different 
from Galactic X-ray binaries, where corona estimates 
range over ^20,000-700,000 km for high i nclination sys- 
tems ([Church fe Balucihska-Church| [2004). while Chan- 
dra grating sp ectroscopy of Cygn us X-2 suggests a corona 
~100,000 km (|Shulz et all 120091) . 

The second type results in optically thin coronae with 
T < kT in and is also well described by a disk blackbody 
+ power law model for som e sources; thes e spectra rep- 
resent an extended corona (|Barnardll2010D . and may be 
related to the steep power law state observed in canoni- 
cal b lack hole binaries (see e.g. iRemillard fc M cClintock 
2006). The electron temperature determines whether or 
not the Comptonized component is well described by a 
power law; if the electron temperature is too close to (or 
within) the observed energy band, then the power law 
fails and more specialized models of Comptonization are 
required (e.g. COMPTT). These coronae are optically 
thin, meaning that we are able to see a portion of the 
disk emission even though when the disk is extended. 
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Extended coronae allow super-Eddington emission while 
keeping the local accretion rate sub-Eddington. 

When disk-blackbody + power law models are applied, 
the power law component sometimes dominates the low 
energy e mission; this has rega r ded as unphysical by som e 
(see e.g. iRoberts et all l2005t iGoncalves fc Sorial 120061 ). 
This led iSteiner et all ()2009f ) to develop the SIMPL con- 
volution model of Comptonization for XSPEC v. 12 that 
results in a si ngle, Comptonized emission component. 
iBarnardl (|2010Q modeled XMM-Newton emission spectra 
of three ULXs in NGC253 (plus the dynamically con- 
firmed BH+Wolf-Rayet binary IC10 X-l) with two com- 
ponent models and convolution models. The convolution 
model was rejected for NGC253 ULX3 and IC10 X-l (and 
the luminosity for IC10 X-l was well sub-Eddington). 
NGC253 ULX1 and NGC253 ULX2 preferred the two 
component model over the convolution model (with lower 
X 2 for the same number of degrees of freedom) . We in- 
fer from these results that the soft excess is real, and 
favor an extended corona (i.e. To < kT n ). We present a 
possible explanation for the soft excess in the discussion. 

XMMUJ004243. 6+412519 is a transient X-ray source 
in M31, identified for the first time in a 2012 January 
15 XMM-Newton ob servation, hereafter observation XI 
(jHenze et al.l l2012af ) ; its lumino s ity wa s found to be 
~2xl0 38 erg s" 1 . IHenze et all (|2012bt ) observed the 
source again twice more with XMM-Newton (observa- 
tions X2-X3); the highest observed luminosity (in X3) 
was 2.5xl0 39 erg s _1 , assuming a two component emis- 
sion model consisting of a power law with photon index 
T = 2.7±0.1, and a disk blackbody with inner disk tem- 
perature kXl n = 0.79±0.05 keV, suffering line of sight 
absorption equivalent to 6.2±0.5xl0 21 H atom cm~ 2 . 
This luminosity is a factor ^2 higher than t he Eddington 
limit fo r a standard 10 Mq black hole, hence IHenze et all 
(|2012bD labeled it an ultra-luminous X-ray source. Since 
this is the second ULX in M31, we will refer to this source 
a s M31 ULX2 

IHenze et all (j2012d ) followed up the XMM-Newton ob- 
servations of M31 ULX2 with five Swift observations (Sl- 
S5) over 2012 February 19 - March 4 that were consistent 
with constant intensity; they combined the spectra, and 
found the summed spectrum to be well described by ei- 
ther a pure 0.88±0.04 keV disk blackbody or a 0.9±0.3 
keV disk blackbody + power law with T fixed to 2.7. A 
further 2012 May 24 Swift observation (S6) resulted in 
a spectrum that was well described by a 0. 66+0^06 keV 
disk blackbody, abs orbed by 4.01^1 x 10 21 atom cm~ 2 
(IHenze et al.ll2012dr>. 

iMiddleton et al I t 20 121 ) discovered strong variable ra- 
dio emission associated with the X-ray source during out- 
burst, and found it likely that the radio source was a jet 
powered by near-Eddington accretion onto a stellar mass 
black hole. 

They also carefully modeled the X-ray spectra from 
X1-X3, and found that an emission model consisting 
of a disk blackbody (diskbb) + Comptonization model 
(comptt) was preferred for X2 and X3 over the disk 
blackbody + power law model used by Henze et al. 
(2012b); however, they do not discuss this in detail. 
For their diskbb + comptt models, kT in ~0.5 keV, T 
~1.2 keV, and the optical depth of the corona was ~11; 
they did not give the electron temperature. The result- 



ing 0.3-10 keV luminosities were considerably lower than 
those obtained by Henze et al. (2012b): 9.8±0.2xl0 38 
and 1.26±0.02x l 39 er g s" 1 for X2 and X3 respectively. 
IMiddleton et al.l (|2012t ) present spectra for M31 ULX2 at 
variou s stages of its spectral evolution. 

The IMiddleton et al.l (pOll) emission mod el represents 
the Tp > kTi n family of models, while the IHenze et al.1 
(|2012b[ ) model represents the To < kTi n model family. 
We will therefore consider the luminosity solutions for 
X2 and X3 from Henze et al. (2012b), as well as those 
of Middleton et al. (2012). 

Two types of X-ray binaries (XBs) exhibit transient be- 
havior. Low mass X-ray binaries may be transient X-ray 
sources due to instabilities in their accretion disks; the 
disk has two stable phases (hot and cold), and an un- 
stable intermediate phase — matter accumulates in the 
disk in the cold phase, and is rapidly du mped onto th e 
compact object in the hot phase (see e.g. Lasota 2001). 
However, the X-rays produced by accretion from the hot 
disk prevent the disk from cooling; the X-ray luminos- 
ity decays exponentially if the whole disk is ionized, and 
linea rly if only part of the disk is ionize d (|King fe R ittcr 
119981) . Meanwhile high mass X-ray binaries with large 
eccentricities may be transient if mass transfer is only 
possible near periastron, or if Be class donor sta rs expe- 
rience ejection events (see e.g. lStella et alJ[l986| ). 

We have been monitoring the central region of M31 for 
the last ^13 years with Chandra, averaging ~1 observa- 
tion per month, looking for X-ray transients. Promising 
examples are followed up with two HST ACS observa- 
tions, the first is taken a few weeks after outburst, and 
the second observation is taken ~6 months later; this al- 
lows us t o identify the counter part via difference imaging 
(s ee e.g.|Barnard et al.ll2012bL and references within). 

Ivan P aradiis fc McClintockl (|1994l ) found an empirical 
relation between the ratio of X-ray and optical luminosi- 
ties of Galactic X-ray binaries and their orbital periods, 
suggestive that the optical emission is dominated by re- 
processed X-rays in the disk; this relation appears to 
be independent of inclination. Their chosen X-ray band 
was 2-10 keV. For an irradiated accretion disk with ra- 
dius a, X-ray luminosity Lx, optical luminosity L opt , and 
temperature T, T 4 cx Lx/a 2 , while the surface bright- 
ness of the disk, S, cx T 2 for t ypical X-ray binaries 
van Paradiis fc McClintockl Il99l . Since L opt cx S.a , 



T T V2 

-^opt CX lv x O, 



2/3 

also a cx P or ' b , where P or b is the orbital 



p eriod. 

van Paradiis fe McClintockl ((1994) defined £ = 
(Tx/T E dd) 1/2 (P OT b/lhr) 2/3 , choosing L EDD = 2.5xl0 38 
erg s" 1 as a normalizing constant, and found 



1.57(±0.24) - 2.27(±0.32) logE. 



(1) 



However, Ivan Paradijs fc McClintockl (|1994[ ) sampled a 
mixture of neutron star and black hole binaries, in var- 
ious spectral states. A cleaner sample was obtained by 
A. Moss et al. (2013, in prep), who used only black hole 
transients at the peaks of their outburst, and found 



M v = 0.84(±0.30) - 2.36(±0.30) logE. 



(2) 



We note that these two relations only differ significantly 
in normalization, caused by black hole X-ray binaries 
having larger disks than neutron star binaries with the 
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same period. We have period estimat es for 12 M31 
transie nts observed by Chandra and HST (Barnard et al. 
I2012bl) . 

We used our first Chandra observation of M31 ULX2 
(2012 February 19, CI) to refine its X-ray position 
to 00:42:43.683 +41:25:18.53 in J2000 co-ordinates, 
with la uncertainties of 0.20" in RA and 0.14" in 
Dec (|Barnard et al.l l2012al ). We then triggered two 
HST/ ACS observations, using the F435W filter that ap- 
proximates the B band. 

In this work we present analysis of our Chandra and 
HST observations of M31 ULX2 (C1-C9 and H1-H2 re- 
spectively), as well as a serendipitous observation in a 
120 ks XMM-Newton observation (PI R. Barnard); we 
refer to this observation as X4. We also re-analyzed the 
Swift observations S1-S6, in order to obtain uncertainties 
on their fluxes; X1-X3 were not available to the public. 
We also give a range of orbital periods estimated for the 
system. 

2. OBSERVATIONS AND DATA REDUCTION 

We observed M31 ULX2 nine times with Chandra 
ACIS over 2012 February - September (C1-C9); unfor- 
tunately, M31 was unobservable for much of March and 
April, resulting in a ^100 day delay between CI and C2. 
Our two HST observations were made on 2012 March 7 
(HI) and September 7 (H2). We scrcndipitously observed 
M31 ULX2 with XMM-Newton for 120 ks on 2012 June 
26 (X4, PI R. Barnard). 

2.1. Optical analysis 

All optical analysis was performed with pc-iraf Revi- 
sion 2.14.1, except where noted. Each HST observation 
included four flat-fielded (FLT) images, and one drizzled 
(DRZ) image. The flat fielded images are corrected for 
instrumental effects, but not background subtracted; the 
total number of counts in each pixel is given. The native 
ACS resolution is com parable to the FWHM of the PSF 
ijFruchter et al.ll2009f ). The drizzled image combines the 
flat-fielded images, removes any cosmic rays, and sub- 
tracts the sky background; it is normalized to give the 
number of counts per second per pixel. We used the DRZ 
images from HI and H2 to create a difference image; 
however, we mainly used the FLT images for our aper- 
ture photometry because the software used (daophot) 
prefers images that include the sky background and mea- 
sure brightness in counts. 

2.1.1. Creating a difference image 

We reprojected the H2 DRZ image into the coordinates 
of the HI DRZ image, to produce an accurate difference 
image. To do this, we first registered the H2 image to 
the HI image with CCMAP, using unsaturated stars that 
were close to the target; this maximized the accuracy of 
the registration at the position of the target. Then, we 
used the iraf task wregister to make the pixel orien- 
tation of the H2 image match that of the HI image. We 
registered the H2 image to the HI image before mapping 
to reduce the noise during image subtraction. The dif- 
ference image was produced by subtracting H2 from HI 
using the ftools task farith. 

2.1.2. Measuring the optical counterpart 



For our first HST observation (HI) we used the 
daophot package released with iraf to obtain the net 
source counts in the FLT images, for a total of Ctot 
counts over T seconds. In particular, we used DAOEDIT 
to obtain the number of source counts from radial profile 
fitting. We converted this to Vega B magnitude via 

B ~ —2.5 log [C't t/T] + ZP, (3) 

having obtained the zero point (ZP — 25.77) from the 
ACS Z ero Point calculatoiQ; we see from iSirianni et al.l 
(|2005l Equation 12 and Table 18) that the conversion 
from F435W counts to B magnitude is within 3er of our 
ZP for B-V = — 0.09 (assuming a typical disk spectrum 
ILiu et all 12001 . 

We can convert from B magnitude to My via 

My = B + 0.09- A H x (1 + 1/3) /1.8 x 10 21 - 24.47, (4) 

where Nh is the line of sight absorption; this accounts 
for the difference in B and V magnitudes of a typical 
accretion disc, a relationship between B band extinction 
and measured line-of-sight absor ption towards the ob- 
ject, and the distance to M31 (see lBarnard et aLll2012bl 
and references within). 

For the second observation (H2) we examined each 
FLT image and found the total number of background 
counts in a circular region with 3 pixel radius at the po- 
sition of the counterpart; this corresponds to ^3 times 
the FWHM width of the PSF. We obtained the 4ct B 
magnitude upper limit from 4 (Ctot) /T. 

2.2. X-ray Analysis 

We extracted spectra from our X-ray observations of 
M31 ULX2 using the appropriate mission-specific soft- 
ware suites: CIAO v.4.4 for C1-C9, HEASOFT v6.13 
for S1-S6, and XMM-Newton SAS v.12.0.1 for X4. We 
used the CALDB version distributed with CIAO. Spec- 
tral analysis was performed with XSPEC vl2.7.1b. 

2.2.1. Chandra analysis 

For each of our Chandra ACIS observations (C1-C9), 
we extracted source and background 0.3-7.0 keV spec- 
tra from circular regions with 20" radii; the high off- 
axis angle meant that the PSF was considerably larger 
than for an on-axis point source. We then created a re- 
sponse matrix using mkacisrmf, and obtained an ancil- 
lary response file from MKARF. The source spectra were 
grouped to give a minimum of 20 counts per bin using 

GRPPHA. 

2.2.2. Swift analysis 

Following the recommended procedure, we extracted 
source spectra from a circular region with 20 pixel (47" ) 
radius, using xselect v2.4b. The background spectra 
were accumulated from a nearby circular region with 40 
pixel (94") radius. We then produced ancillary response 
files using XRTMKARF, using the included exposure map 
to account for variations in effective area over the obser- 
vation. Swift response files are not generated on the fly; 
instead, the correct response is found by using the QUZ- 
CIF tool, and copied into the working directory. Spectral 

1 http:/ /www.stsci.edu/hst/acs/analysis/zeropoints/zpt.py 
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Figure 1. Details of two HST/ ACS observations of M31 ULX2 using the F435W filter, along with the difference image. North is up, east 
is left, a 1" scale is indicated. The 3c uncertainty in the X-ray position is indicated by an ellipse. A white object in the difference image 
is brighter in the first observation than in the second one; the counterpart to M31 ULX2 is easily identified. 



channels 0-29 were labeled bad, and the spectra were 
grouped to give at least 20 counts per bin using GRP- 

PHA. 

2.2.3. XMM-Newton analysis 

Due to the high intensity of M31 ULX2, we only- 
used data from the pn instrument. XMM-Newton 
observations often experience intervals of greatly in- 
creased background levels. We searched for such flar- 
ing intervals by creating a lightcurve with xmmse- 
lect, using the expression "(PATTERN==0)&&(PI in 
[10000:12000])&&(FLAG==0)" and 100 s bins; we fil- 
tered out all intervals where the pn rate was >0.4 count 
s -1 using tabgtigen. Flaring was substantial, with 75 
ks of good time out of a 120 ks total. 

We extracted the source spectrum from a region that 
was optimized by the analysis software; this was a cir- 
cle with radius ~24" . A circular background region was 
chosen to be near the source, on the same chip, and at a 
similar off-axis angle; its radius was ~40". These spec- 
tra were filtered by the good time interval, and by the 
expression "(PATTERN<=4)&&(FLAG==0)" . We ob- 
tained a corresponding response file with rmfgen and 
an ancillary response file with mkarf. The source spec- 
trum was grouped to a minimum of 20 counts per bin, 
and energies outside the range 0.3-10 keV were excluded. 

2.3. Locating the X-ray source 

We were unable to register the Chandra position of 
M31 ULX2 directly with the HST observations due to a 
lack of known X-ray sources in the HST field. Instead we 
used 27 X-ray bright globular clusters (GCs) to register 
a combined ~350 ks ACIS image (supplied by Z. Li) to 
the B band Field 5 image of M31 provided by the Local 
Galaxy group Survey (LGS) (jMassev et all 12006). We 
used PC-IRAF v2.14.1 to perform the regis tration, follow- 
ing the same procedure as described in Barnard et al. 
(|2012bf ). The X-ray and LGS positions of the 27 X- 
ray bright GCs were determined using imcentroid; the 
equivalent FK5 coordinates were calculated for the X- 
ray and optical position of each GC using xy2sky v2.0, 
distributed with ftools. The X-ray positions of each 
GC were altered to match the LGS positions, allowing 
the registration of the merged Chandra image to Field 5 
using the iraf task CCMAP. This registration yielded la 
position uncert ainties of 0.11" in R.A., and 0.09" in Dec 
ijBarnard et"alll2012bD . 

We then registered our first HST observation (HI) to 
the Field 5 observation in the same way using bright, 
unsaturated stars. Similarly, we registered our brightest 



Chandra observation of M31 ULX2 (CI) to the merged 
Chandra image using a selection of bright X-ray sources. 
The final uncertainties in the X-ray position of M31 
ULX2 combine the position uncertainties in the X-ray 
image, the uncertainties in registering Observation CI 
to the merged Chandra image, and the uncertainties in 
registering the merged Chandra image to the M31 Field 
5 LGS image. The uncertainties in registering the HST 
observation were negligible. 

3. RESULTS 

3.1. Optical analysis 

3.1.1. Difference imaging 

The difference image produced from the HI and H2 
images revealed an optical counterpart at 00:42:43.699 
+41:25:18.54 (J2000), with la uncertainties of 0.21" RA 
and 0.18" Dec with respect to the LGS M31 Field 5 
image. We present details of each HST observation in 
Fig. [H along with the difference image. 

3.1.2. Photometry 

We first obtained source counts in HI for the ULX2 
counterpart from FLT images 2 and 3 using PSF fitting 
via DAOEDIT: 1 and 4 were contaminated. For Frame 2, 
we obtained 1010 net source counts over 1255 s; the sky 
mean was 141 count pixel -1 , with a = 33 count pixel -1 . 
For Frame 3, we obtained 1198 net source counts over 
1410 s (sky mean = 166 count pixel -1 , a = 39 count 
pixel -1 ). Individually, these frames gave a B magni- 
tude of 26.03±g;g| and 25.95^;^ respectively. Combin- 
ing these frames yielded B = 25.97^q'q2- For consistency, 
we also checked the B magnitude derived from the DRZ 
image; this yielded a source intensity of 0.867 count s -1 
over 5450 s (sky mean = —0.026 count s -1 pixel -1 , a 
= 0.011 count s -1 pixel -1 ); the DRZ image gave B = 
25.924Jq 1 5 17 , consistent with the magnitude from our 
FLT images. We will use our B magnitude from the 
combined FLT frames (B = 25-97±g;g|) for the rest of 
this work. 

We saw no sign of the counterpart to M31 ULX2 in 
observation H2. Instead we estimated the Aa upper limit 
to the B magnitude from the background counts; we note 
that 3 pixel radius circles used may be contaminated by 
neighboring stars, but this is not too important since 
we are only interested in upper limits. FLT images 1-3 
yielded a total of 6829 counts over 3861 seconds; image 
4 was contaminated by a cosmic ray. Hence, B > 28.4 at 
the Aa level. 
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The optical luminosity of the M31 ULX2 counterpart 
varied by at least ~1 order of magnitude between ob- 
servations. Using the value for Ah obtained from X4, 
My = -0.91 for HI, and >+1.5 for H2; the counter- 
part was more luminous (in terms of absolute magni- 
tude) than any other transi ent counterpart found in M31 
to date, apart from ULX1 (Bar nard et al.ll2012bD ; how- 
ever, it was also fainter (in apparent magnitude) than 
any other counterpart due to the high absorption (see 
below). The known counterparts of HMXBs in the Small 
Magellanic Cloud have absolute values —6 < My < — 1; 
for the known BH HMXBs (Cygnu s X-l, LMC X-l and 
LMC X-3), -6.5 < M v < -1.5 (|Barnard et all l20ilL 
and references within). All of these high mass systems 
are optically brighter than the M31 ULX2 counterpart 
even in outburst; hence a low- mass donor is likely. 

3.2. X-ray analysis 

We summarize the X-ray properties of M31 ULX2 dur- 
ing all the observations that we analyzed in Table [TJ 
For each observation we give the time since the first 
observation, the absorption, as well as the disk black- 
body temperature and/or power law photon index de- 
pending on the spectral state; the hard state was repre- 
sented by a power law (xspec model WABS*CFLUX*POW) 
with photon index 1.7 (as is typical), the thermally dom- 
inated state was represented by WABS*CFLUX*DISKBB, 
while the steep power law state was represented by 
WABS*(cflux*diskbb + CFLUX*pow). Finally we give 
the x 2 /dof, luminosity and spec tral state; luminosi- 
ties a ssume a distance of 780 kpc (jStanek fc Garnavichl 
1998). Values quoted without uncertainties are frozen. 
All uncertainties in results presented in this work are 
quoted at the la level. 



Unfolded Spectrum 




Energy (keV) 

Figure 2. The 0.5-10 keV S4 XRT spectrum fitted with the best 
Scenario B model. The disk blackbody and power law components 
are represented by dashes and dots respectively. A color version is 
provided in the electronic edition. 

The spectra of X2 and X3 have been interpreted in 
two ways. Henze et al. (2012b) used a disk blackbody 
+ power law emission model, and obtained a 0.3-10 
keV luminosity of ~2.5xl0 39 erg s _1 for X3. Alterna- 
tively, Middleton et al. (2012) used a disk blackbody + 
compTT emission model with a 0.3-10 keV luminosity 
for X3 of 1.26±0.02xl0 39 erg s _1 . We note that the fit 
obtained by Middleton et al. (2012) suggests a compact, 



optically thick, opaque corona, while the Henze et al. 
mo del represents an ext ended, optically thin corona (see 
e.g. lBarnard et alj[201lh : we shall label these Scenario A 
and Scenario B respectively. 

3.2.1. Re- analyzing the Swift data 

Following Henze et al. (2012c), we fitted S1-S5 simul- 
taneously, with all of the parameters free to vary but 
tied between observations; we obtained similar results. 
However, we note that the number of degrees of free- 
dom in our fits were substantially smaller than for the 
fits quoted by Henze et al. (2012c); this is probably due 
to different binning strategies. Our best simultaneous 
fit for a disk blackbody emission model yielded Ah = 
2.87±0.15 x 10 21 atom cm" 2 , and kT in = 0.877±0.016 
keV; x 2 / dof = 377/332. When fitting a disk blackbody 
+ power law model, the photon index was frozen at 2.7, 
following lHenze et al.l (|2012b| ): A H = 3.9t^ x 10 21 atom 
cm- 2 , kTi n = 0.85±0.02 keV, and x 2 /dof = 374/331. 

The null hypothesis probabilities for the model fits to 
S1-S5 are 9% for Scenario A, and 10% for Scenario B. We 
caution that there is no statistical requirement for the ad- 
dition of the power law component to these relatively low 
quality spectra; our motivation for including Scenario B 
is the discrepancy between the Scenario A lightcurve and 
our optical results, c.f. the known behavior of Galactic 
transients (see below) . We also note that our unabsorbed 
luminosities for Scenario B assume a particular power 
law index (following Henze et al., 2012bc); Scenario B 
luminosities may therefore be subject to systematic un- 
certainties; however, our X4 spectrum is consistent with 
T = 2.7 at the la level. 

The 0.3-10 keV luminosity for the disk blackbody 
model for S1-S5 was l.ll±0.02x 10 39 erg s -1 , while 
the two-component model yielded a luminosity of 
1.5±0.3xl0 39 erg s . These two luminosities are con- 
sistent, due to the large uncertainties in the luminosity 
from the two-component model. In physical terms, the 
extended corona would appear brighter because cool pho- 
tons in the outer disk are up-scattered into the observed 
energy range. 

In Fig. [2] we present the best fit unfolded Scenario B 
spectrum multiplied by energy for S4 (the highest qual- 
ity Swift spectrum); the disk blackbody and power law 
components are represented by dashes and dots respec- 
tively. We see that the power law component dominates 
below ~1 keV; this is physically possible for an extended 
corona but not for a compact corona. The best fit x 2 /dof 
= 98/87. Setting the power law contribution to zero in- 
creases the % 2 /dof to 148/88, demonstrating that the 
power law contribution is significant in this model; how- 
ever, a single disk blackbody emission model fits this 
spectrum just as well, with x 2 /dof = 99/88. 

Our analy s is of S6 supports that analysis reported by 
I Henze et al.l (|2012dl ). Fitting an absorbed disk black- 
body model yielded A H = 3.5±0.9xl0 21 atom cm -2 , and 
kT in = 0.62±0.05 keV; x 2 /dof = 19/18 with a 0.3-10 keV 
luminosity of 4.8±0.6xl0 38 erg s _1 . 

3.2.2. The X-ray evolution of MSI ULX2 

Both scenarios agree on the following. The X-ray emis- 
sion of M31 ULX2 evolved considerably during the out- 
burst: it appears to have evolved from the hard state 
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Table 1 

Summary of our spectral fitting results for M31 ULX2. For each observation we give the time from first observation, absorption 
normalized to 10 21 atom cm -2 , disk blackbody temperature, and power law photon index; we then give the corresponding x 2 /dof, 0.3—10 
keV luminosity normalized to 10 37 erg s _1 , and spectral state. For CI and S1-S5 we give fits pertaining to Scenario A (TD) and B 
(SPL); S1-S5 were fitted simultaneously with all parameters tied. The quality of the spectrum is indicated by the degrees of freedom 
(dof), since each spectral bin contains at least 20 counts. For C8-C9, we assume a hard state emission model, with JVh = 3.37xl0 21 atom 

cm~ 2 , and T = 1.7. All uncertainties are quoted at the ler level. 



Obs 


T-T 




kT in / keV 


r 


X^/dof 


£37 


State 


SI 


34.4 


2.87±0.15 


0.877±0 


016 


— 


377/332 


111±2 


TD (A) 


CI 


34.6 


2.4±0.5 


0.97±0 


05 


— 


131/127 


106±3 


TD (A) 


S2 


39.2 


2.87±0.15 


0.877±0 


016 


— 


377/332 


111±2 


TD (A) 


S3 


46.6 


2.87±0.15 


0.877±0 


016 


— 


377/332 


111±2 


TD (A) 


SI 


47.4 


2.87±0.15 


0.877±0 


016 


— 


377/332 


111±2 


TD (A) 


S5 


48.4 


2.87±0.15 


0.877±0 


016 


— 


377/332 


111±2 


TD (A) 


SI 


34.4 




0.85±0 


02 


2.7 


374/331 


150±30 


SPL (B) 


CI 


34.6 


4-5+^ 


0.94±0 


04 


2.7 


127/126 


175±50 


SPL (B) 


S2 


39.2 


3 9 + °i 

° -i7 — 0.5 


0.85±0 


02 


2.7 


374/331 


150±30 


SPL (B) 


S3 


46.6 


3 9+2' 4 


0.85±0 


02 


2.7 


374/331 


150±30 


SPL (B) 


S4 


47.4 


n n + 0.4 
°' a -0.5 


0.85±0 


02 


2.7 


374/331 


150±30 


SPL (B) 


S5 


48.4 


3 n+0.4 


0.85±0 


02 


2.7 


374/331 


150±30 


SPL (B) 


HI 


50.9 
















S6 


130.0 


2.9±0.7 


0.64±0 


05 




20/19 


45±4 


TD 


C2 


131.9 


5. Oil. 7 


0.62±0 


07 




14/15 


47±8 


TD 


C3 


137.5 


3.18±0.10 


0.656±0 


007 




257/222 


53.0 ±0.7 


TD 


C4 


142.7 


3.24±0.11 


0.645±0 


007 




196/210 


48.7±0.7 


TD 


C5 


148.7 


3.23±0.11 


0.641±0 


007 




200/213 


45.0±0.6 


TD 


C6 


157.8 


3.2±1.2 


0.61+.0 


05 




28/23 


33±5 


TD 


XI 


163.3 


3.37±0.07 


0.577+.0 


005 


n q+0.4 


684/634 


40.7+.1.0 


TD 


C7 


186.3 


3.5±1.4 


0.58±0 


08 




18/19 


20±3 


TD 


C8 


216.0 


3.37 






1.7 


N/A 


2.6±0.4 


HS 


H2 


234.9 
















C9 


240.6 


3.37 






1.7 


N/A 


0.7±0.3 


HS 



(HS) to a two component emission state during the rise 
(X1-X3), either the "ultra-luminous" state or something 
resembling the the steep power law (SPL) state; we note 
that the true peak may not have been observed. It then 
went to the thermally dominated (TD) state in decline 
(X4, C2-C7, S6). A sudden drop in intensity at the end 
of the outburst (C8-C9) suggests a transition from the 
TD state (where the most of the bolometric luminosity 
is in the observed band) to the HS state (where much of 
the power can be emitted at much higher energies). 

The difference between Scenarios A and B is in the 
interpretation of the spectra from CI and S1-S5. These 
spectra are consistent with being in a SPL-like state (disk 
blackbody + power law spectra), or in the TD state (disk 
blackbody emission only). The TD spectra for CI and 
S1-S5 are only consistent with Scenario A; the SPL spec- 
tra are consistent with both scenarios due to large un- 
certainties. 

For Scenario A, the shape of the 0.3-10 keV lightcurve 
is more linear than exponential, losing ^5xl0 36 erg s _1 
per day. However, Scenario B favors exponential decay, 
with a 0.3-10 keV e-folding time of ~80 days. Since 
X-ray transient decay is expected to be exponential or 
linea r depending on whethe r or not the disk is fully ion- 
ized, iShahbaz et. al.l ([1998) studied the dependence of 
lightcurve shape on orbital period and peak X-ray lumi- 
nosity. For black hole transients with peak X-ray lumi- 
nosities >10 39 , they expect exponential decay even for 
systems with periods >300 hr. For Scenario A to be cor- 
rect, the optical counterpart is expected to have My < 
—5 according to either Equation 1 or 2, and B <21 ac- 
cording to Equation 4; this is considerably brighter than 



the observed mag nitude, B = 25.951^, hence we find 
Scenario A unlikely. 
We present a lightcurve of 0.3-10 keV luminosities for 
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Figure 3. Unabsorbed luminosity lightcurve of M31 ULX2 in the 
0.3—10 keV band. Black open squares represent XMM-Newton lu- 
minosities taken from Hcnzc et al (2012ab); no uncertainties were 
given. Gre y open squares r e presen t the luminosities for X1-X3 
obtained by Middlcton et al. (2012). The re-analyzed Swift obser- 
vations are represented by open diamonds. Our Chandra ACIS 
observations are shown as filled circles, and our XMM-Newton ob- 
servation is represented by a filled square; uncertainties are quoted 
at a 1<t level. For CI and SI— S5, gray and black symbols represent 
spectra for Scenarios A (TD) and B (SPL) respectively. HST ob- 
servation times are indicated by downward arrows. Vertical lines 
split the lightcurve by inferred spectral state: hard state (HS), 
thermally dominated (TD), and steep power law (SPL), following 
Rcmillard & McClintock ( 2006). The green and blue lines repre- 
sent linear and exponential respectively for Scenario A, while the 
red line represents exponential decay for Scenario B. 



M31 ULX2 with Chandra, HST and XMM-Newton 
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M31 ULX2 in Fig. |31 created using the results of Henze 
et al. (2012ab) and Middleton et al. (2012) in addition 
to our Chandra, Swift, and XMM-Newton results. We 
present only one datum for observations that both sce- 
narios agree on: C2-C9, X4, and S6. Grey and black 
data correspond to Scenarios A and B respectively for 
X1-X3, CI, and S1-S5. X1-X3 are represented by open 
squares, X4 by a filled square, C1-C9 by circles, and Sl- 
S6 by diamonds. Times of HI and H2 are indicated by 
downward arrows. 

We used our X4 0.3-10 keV pn spectrum of M31 
ULX2 to get the best measurement of the absorption; 
we obtained ~62,000 net source counts. An absorbed 
disk blackbody provided an acceptable fit (x 2 /dof = 
689/636). This was slightly improved by adding a power 
law component (% 2 /dof =684/634): absorption Njj = 
3.371qq5 x 10 21 atom cm~ 2 , inner disk temperature kT ln 

= 0.577±0.005, photon index T = 2.3±£f , quoting la un- 
certainties. The total 0.3-10 keV unabsorbed luminosity 
was 4.07±0.10xl0 38 erg s" 1 , with 3.95±0.04x 10 38 erg 
s _1 contributed b y the disk blackbody; this corres ponds 
to the TD state (jRemillard fc McClintockl [2006). We 
present this best fit to the pn spectrum in Fig. |4] 

While most of the absorption measurements are con- 
sistent with that from X4, we no te that the spectral fits 
obtained by iHenze et al.l (I2012bf) for X2 and X3 sug gest 
significantly higher absorption. IHenze et al.l (|2012bl ) use 
a power law component to represent the Comptonized 
emission component; this could overestimate the num- 
ber of soft photons in the unfolded spectrum, thereby 
increasing the absorption required to produce observed 
spectrum. However, observations S1-S5, and CI were 
also modeled with the same power law component in Sce- 
nario B, and yielded Nh values consiste nt with X4. Fur- 
therm ore, the absorption obtained by iMiddleton et al.1 
(2012) for X2 is 8a higher than our absorption for X4. 
Hence the increase in absorption during outburst is likely 
to be real; this extra absorbing material may have been 
ejected during the outburst, and dissipated by the time 
of X4. 

We also note that the absorption measured in the si- 
multaneous fitting of S1-S5 in Scenario A is 3. Oct lower 
than that measured in X4, further evidence against Sce- 
nario A. 

3.3. Estimating the orbital period 

3.3.1. Scenario A 

In Scenario A, the 0.3-10 keV luminosity of M31 ULX2 
at the time of HI was ^9.7xl0 38 erg s _1 , assuming a 
decrease of ~5xl0 36 erg s _1 day -1 . Assuming the CI 
TD spectrum (Table 1), this corresponds to a 2.0-10 keV 
luminosity — 4.5xl0 38 erg s -1 . Since M31 ULX2 is in 
the TD state for Scenario A, we may use Equation 2 to 
estimate the orbital period. The suggested orbital period 
is 11±8 hr at the la level, including uncertainties in the 
relation; this is clearly shorter than the >300 hr expected 
from the linear decay ()Shahbaz et. al.|["l998f) . 

3.3.2. Scenario B 

For Scenario B we estimated the 2-10 keV luminosity 
of M31 ULX2 at the time of HI to be -4.4 xlO 38 erg s~ f . 
Using the relation specific to black holes, we obtained an 



orbital period of 8.5±0.4 hr, not including uncertainties 
in the relation. 

However, we note that M31 ULX2 is likely to have 
been in the steep power law state during HI, whereas 
the black hole period relation (Equation 2) is calibrated 
for thermally dominated spectra; the X-ray to optical 
ratio is likely to be different for different spectral states 
in the same system. Using Equation 1 results in an or- 
bital period of 28.5±1.3 hr; given that the sample used 
in deriving Equation 1 was dominated by neutron stars, 
and that black hole binaries with the same size disk have 
shorter periods than neutron star binaries, we consider 
this period to be some sort of upper limit. 

A period for M31 ULX2 in the range ~9-30 hr is con- 
sistent with the period estimated for the first ULX in 
M31 (T9 in Barnard et al., 2012b, P = 18^ hr). 

4. DISCUSSION AND CONCLUSIONS 

XMMUJ004243.6+412519 (M31 ULX2) is the second 
ULX to be observed in M31; both of these systems are 
transient. We have combined analysis of our Chandra, 
HST, and XMM-Newton observations of M31 ULX2 with 
the reanalysis of Swift data and reported results of Henze 
et al. (2012abcd) and Middleton et al. (2012) to build a 
picture for the system. 

We have identified an optical counterpart by its change 
in B magnitude from 25.97 to >28.4 (My from -0.91 to 
>+1.5). Hence the optical emission is dominated by re- 
processed X-rays from the disk. Since even the peak opti- 
cal emission is fainter than known HMXB counterparts, 
a low mass donor is expected. We estimate an orbital 
period of ~9-30 hr; this period is comparable with that 
estimated for the first M31 ULX. 

Two solutions were presented for fitting the spectra 
from the highest luminosity observations. Middleton et 
al. (2012) used a disk blackbody + Comptonization 
model that represented a typical "ULX state", with a 
cold, optically th ick corona that obscu res the inner re- 
gions of the disk (|Gladstone et al.ir2009D ; we refer to this 
as Scenario A. Henze et al. (2012b) used a disk black- 
body + power law model that represents a hot, optically 
thin corona tha t extends over much of the accretion disk 
(Barnard 2010); we refer to this as Scenario B. Scenario 

Unfolded Spectrum 




0.5 1 2 

Energy (keV) 



Figure 4. Unfolded 0.3-10 keV XMM-Newton pn spectrum of 
M31 ULX2 from 2012 June 26, fitted with a disk blackbody (red 
dotted line) and a power law component (blue dashed line). The 
disk blackbody component contributes ~97% of the total 0.3—10 
keV emission. A color version is provided in the electronic edition. 
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A results in a linear decay lightcurve, while Scenario B 
yields an exponential decay lightcurve. For X-ray lu- 
minosities >10 39 erg s _1 , linear decay is only expecte d 
for systems with periods 3>300 hr ([Shahbaz et. al.lll998h . 
Since our orbital period estimate is ~9-30 hr, we favor 
Scenario B for this system. However, we caution that 
there is no statistical evidence in the X-ray spectra that 
Scenario B is superior. Furthermore, the Scenario B lu- 
minosities assume T=2.7, so the luminosities may be sub- 
ject to systematic uncertainties. 

Since this work provides further support for spectral 
models where the Comptonized compo nent dominates 
the low energy emi ssion (on top of e.g. Barnard] 120101 : 
iBarnard et atll2l)l~ll ) . we consider how this could be pos- 
sible. We are assuming unsaturated, inverse Compton 
scattering of cool photons on hot electrons in the corona; 
the mean number of scatterings experienced by the pho- 
tons is much less than the number of scatterings required 
to bring the photons up to the electron energy (see e.g. 
iSunvaev fc TitarchukJ Il980l and references within). If 
the soft excess is due to Comptonization of disk photons, 
the lower energy photons must be scattered more often 
than the high energy photons. This could be achieved if 
the corona density decreased in the inner disk region due 
to increased radiation pressure; the high energy photons 
in the inner disk would experience fewer scatterings than 
the low energy photons in the outer disk. 

The spectral evolution of M31 ULX2 was unusual in 
that it went from something like the SPL state to the TD 
state during the decay. However, we note t hat the spec- 
trum re ported for X3 (the observed peak) bv lHenze et al.l 
(|2012bD is similar to the spectra exhibited by NGC253 
ULX1, which in turn re sembles the bl ack hole + Wolf- 
Rayet binary IC10 X-l (|Barnardll20l"ol ). The properties 
of those systems suggested extended coronae, and we 
proposed that many ULXs could be explained by stel- 
lar mass black hole systems with truly super-Eddington 
accretion overall, while the local accretion rate remained 
sub-Eddington. It is intriguing that in all of the observa- 
tions where M31 ULX2 is super-Eddington, its spectra 
were reminiscent of the SPL state. However M31 ULX2 
does appear to be consistent with a low mass X-ray bi- 
nary containing a ~10 M© black hole. 

It is interesting to compare the outbursts of the two 
transient ULXs in M31. The outburst of the first tran- 
sient ULX appears to have been more traditional, going 
from the HS to the TD state during the rise, then to the 
SPL and back to the HS state during decay; indeed it 
even exhibited a TD state at a 0.3-10 keV lum inosity 
of 1.75±0.07xl0 39 erg s" 1 (INooraee et alJ[20H . M31 
ULX1 is also a low mass X-ray binary, but the mass of 
the black hole appears to be somewhat higher than usual. 
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